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The production of D mesons in the scattering of transversely polarized protons off unpolarized 
protons at RHIC offers a clear opportunity to gain information on the Sivers gluon distribution 
function. D production at intermediate rapidity values is dominated by the elementary gg cc 
channel; contributions from qq — > cc s-channel become important only at very large values of a;^. 
In both processes there is no single spin transfer, so that the final c or c quarks are not polarized. 
Therefore, any transverse single spin asymmetry observed for D's produced in p^p interactions 
cannot originate from the Collins fragmentation mechanism, but only from the Sivers effect in the 
distribution functions. In particular, any sizeable spin asymmetry measured in p^p —* D X at 
mid-rapidity values will be a direct indication of a non zero Sivers gluon distribution function. We 
study the p^p — > _D X process including intrinsic transverse motion in the parton distribution and 
fragmentation functions and in the elementary dynamics, and show how results from RHIC could 
allow a measurement of fg^^f . 

PACS numbers: 13.88.-l-o, 12.38.Bx, 13.85.Ni 



I. INTRODUCTION AND FORMALISM 



Within the QCD factorization scheme, the cross section for an inclusive large pt scattering process between hadrons, 
like pp —^ h + X, is calculated by convoluting the elementary partonic cross sections with the parton distribution 
functions (pdf's) and fragmentation functions (ff's). These objects account for the soft non-perturbative part of 
the scattering process, by giving the probability density of finding partons inside the hadrons (or hadrons inside 
fragmenting partons) carrying a specific fraction x (or z) of the parent light-cone momentum. The parton intrinsic 
motion - demanded by uncertainty principle and gluon emission - is usually integrated out in the high energy 
factorization scheme, and only partonic coUinear configurations are considered. However, it is well known that the 
quark and gluon intrinsic transverse momenta k± have to be taken into account to improve agreement with data on 
unpolarized cross sections at intermediate energies Moreover, without intrinsic k± one would never be able to 

explain single spin asymmetries within QCD factorization scheme; several large single spin asymmetries have been 
observed 0,0, Q , which in a coUincar configuration are predicted to be either zero or negligibly small. Although not 
rigorously proven in general |^ , the usual factorized structure of the coUinear scheme has been generalized with 
inclusion of intrinsic fcj^, so that the cross section for a generic process A B ^ C X reads: 

dcr = ^ fa/A{xa,k±a) fb/B{xb,k±b)<^da''^^''-{xa,xt,k±a,k±b)<^Dc/c{z,k±c) ■ (1) 

a.b,c 

The pdf's and the ff's are phenomenological quantities which have to be obtained - at least at some scale - from 
experimental observation and cannot be theoretically predicted. The pdf's of unpolarized nucleons, q{x) = fq/p{x) = 
fi{x), are now remarkably well known; one measures them in inclusive deep inelastic scattering processes at some 
scale, and, thanks to their universality and known QCD evolution, can use them in different processes and at different 
energies. The k± dependence of q{x, k±) is usually assumed to be of a gaussian form, and the average k± value can 
be fixed so that it agrees with experimental data. Notice that, in our notations, a hat over a pdf or a ff signals its 
dependence on fcj^ and k± = \k±\. 

When considering polarized nucleons the number of pdf's involved grows and dedicated polarized experiments 
have to be performed in order to isolate and measure these functions. We have by now good data on the pdf's of 
longitudinally polarized protons - the helicity distribution Aq(x) = gKx) - but nothing is experimentally known on 
the transverse spin distribution - the transversity function ATq{x) = Sq{x) = h\{x). The situation gets much more 
intricate when parton intrinsic transverse momenta are taken into account. Many more distribution and fragmentation 
fmictions arise, like the Sivers function /S.^ f{x,k±) (x f^{x,k±) which describes the probability density of 

finding unpolarized partons inside a transversely polarized proton; similarly, the Collins fragmentation function Q 



2 



gives the number density of unpolarized hadrons emerging in the fragmentation of a transversel y p olarized quark. 
These are the functions which could explain single spin asymmetries in terms of parton dynamics |lOl lll| . 

One of the difficulties in gathering experimental information on these new spin and k± dependent pdf's and ff's 
is that most often two or more of them contribute to the same physical observable, making it impossible to estimate 
each single one separately. 

In Ref. |12| | it was shown how properly defined single spin asymmetries in Drell-Yan processes depend only on the 
Sivers distribution function A'^f{x,k±) of quarks (apart from the usual known unpolarized quark distributions). In 
Ref. 13] it has been suggested to look at back-to-back correlations in azimuthal angles of jets produced in p^p RHIC 
interactions in order to access the gluon Sivers function. We consider here another case which, again, would isolate 
the gluon Sivers effect, making it possible to reach direct independent information on fg/pr{x,k±). 

Let us consider the usual single spin asymmetry 
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(2) 



for p^p DX processes at RHIC energy, ^/s = 200 GeV. These D mesons originate from c or c quarks, which at 
LO can be created either via a qq annihilation, qq — )■ cc, or via a gluon fusion process, gg cc. The elementary 
cross section for the fusion process includes contributions from s, t and it-channels, and turns out to be much larger 
than the qq — > cc cross section, which receives contribution from the s-channel alone. Therefore, the gluon fusion 
dominates the whole p^p — > DX process up to xp ~ 0.6. Beyond this the qq — > cc contribution to the total cross 
section becomes slightly larger than the gg — > cc contribution, due to the much smaller values, at large x, of the gluon 
pdf, as compared to the quark ones (see Fig. 1). 

As the gluons cannot carry any transverse spin the elementary process gg — *■ cc results in unpolarized final quarks. 
In the qq —^ cc process one of the initial partons (that inside the transversely polarized proton) can be polarized; 
however, there is no single spin transfer in this s-channel interaction so that the final c and c are a gain not polarized. 
One might invoke the possibility that also the quark inside the unpolarized proton is polarized jl4| . so that both 
initial q and q are polarized: even in this case the s-channel annihilation does not create a polarized final c or c. 
Consequently, the charmed quarks fragmenting into the observed D mesons cannot be polarized, and there cannot be 
any Collins fragmentation effect [see further comments after Eq. ©]. 

Therefore, transverse single spin asymmetries in p^p DX can only be generated by the Sivers mechanism, namely 
a spin-fcx asymmetry in the distribution of the unpolarized quarks and gluons inside the polarized proton, coupled 
respectively to the unpolarized interaction process qq — > cc and gg — *■ cc, and the unpolarized fragmentation function 
of either the c or the c quark into the final observed D meson. That is ^hj: 



EodaP^P^^^ EodaP^P^^^ 



+ 



d^PD d^Pn 
dxa dxb dz d^k±a d^k±b d^ko 5{kD ■ Pc) S{s + i + u — 'i-mq) C{xa,Xb, z, k^) 

^ fq/pr{Xa,k±a) fq/p{xb,k±b) {Xa, Xb, k±a, k±b, kn) Do/Qiz^ko) 

^ " da39~*QQ 

^ fg/pr{Xa,k_l_a) fg/p{xb,k±b) {Xa, Xh, k^a, k^bi ko) DD/Q{z,kD) 



(3) 




where q = u,u,d,d, s,s and Q = c or c, according to whether D — _D+, D° or D — D^ , D . Notice that z is the 
light-cone momentum fraction along the fragmenting parton direction, identified by p^, z — P^/pt ■ Throughout the 
paper we choose XZ as the D production plane, with the polarized proton moving along the positive Z-axis and 
the proton polarization "f along the positive y-axis. In such a frame k±a and k±b have only X and Y components, 
while ko has all three components; the function S{kD ■ p^) ensures that the integral over kjj is only performed along 
the appropriate transverse direction, fc^^i, that is the transverse momentum of the produced D with respect to the 
fragmenting quark direction. The factor C contains the flux and relevant Jacobian factors for the usual transformation 
from partonic to observed meson phase space, which, accounting for the transverse motion, reads [Tsl Il6||: 
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Notice that for coUinear and massless particles this factor reduces to the famihar s/nz^. The Sivers distribution 
functions Q for quarks and gluons are defined by 

^ /a/pT {^ai fc_La) — /a/pT (^aj fe±a) ^ /a/pi {^a, fc±a) — /a/pT {^a; fc±a) ~ /a/pT {^a, ~fe±a) j (5) 

where a can either be a Ught quark or a gluon. Similarly, Dd/q{z, k^o) is the probability density for a quark Q to 
fragment into a D meson with light-cone momentum fraction z and intrinsic transverse momentum k^D. 

The heavy quark mass toq is taken into account in the amplitudes of both the partonic processes and the resulting 
elementary cross sections are: 



dl P 9 



(2rf + 2t| + p) , 



di S2 8 y2,TiT2 "J ■ 4riT2j ' 

where ti_2 and p are dimensionless quantities defined in terms of the partonic Mandelstam variables s, t and u as 

rrin — i 



(6) 
(7) 



T2 = 
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The denominator of ^at, Eq. ||2J), is analogously given by 

, . ,1 EodaP^P^^^ EodaP^P^^^ EodaPP^^^ 
rf-^Pu d'^Po d-^Po 

= 2 / da;}, dz d'^k±a d^k±i, d^ko 5{kD ■ Pc) 5{s + t + ii — "^itiq) C{xa, Xb, z, k^) 



(9) 



E 



di 



{xa,Xb, kj_a, kj_b, ko) Du/q{z, ko) 



fg/p{Xa, k±a) fg/p{xb, k±b) {Xa, Xb, k±a,k±b, ko) Dd/q{z, ko) 

In Eqs. and lO we consider intrinsic transverse motions in the distributions of initial light quarks, in the 
elementary process and in the heavy quark fragmentation function, i.e. we consider a fully non planar configuration 
for the partonic scattering. This has two main consequences: on one side, taking into account three intrinsic transverse 
momenta makes both the kinematics and the dynamics highly non trivial; on the other side it generates a large 
number of contributions, other than the Sivers effect, which originate from all possible combinations of fcj^ dependent 
distribution and fragmentation functions, each weighted by a phase factor (given by some combination of sines and 
cosines of the azimuthal angles of the parton and final meson momenta) . This topic deserves a full treatment on its 
own, which will soon be presented in Refs. |l7l Il8l|. At present, we only point out that we have explicitely verified 
that all contributions to the p^p DX single spin asymmetry from k± dependent pdf's and ff's, aside from those 
of the Sivers functions, are multiplied by phase factors which make the integrals over the transverse momenta either 
negligibly small or identically zero. Therefore, they can be safely neglected here. 

Eq. |j3J| shows how depends on the unknown Sivers distribution function; as all other functions contributing to 
Ak[ are reasonably well known (including the fragmentation function Djj/q 19]) a measurement of An should bring 
direct information on fg/pf and, to some extent, A^/^/pT- 



II. NUMERICAL ESTIMATES 



So far, all analyses and fits of the single spin asymmetry data were based on the assumption that the gluon Sivers 
function A^/g/j,? is zero. RHIC data on A^ in p^p DX will enable us to test the validity of this assumption. In 
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FIG. 1: The unpolarized cross section for the process pp — » DX at ^/s = 200 GeV, as a function of Ed and pr at fixed 
pseudo-rapidity rj = 3.8 (a), and as a function of xf at fixed transverse momentum pT = 1.5 GeV/c (b), calculated according 
to Eqs. Q and HIOI I. The solid line is the full cross section, whereas the dashed and dotted lines show the qq — > cc and gg — > cc 
contributions separately. 



fact, as the gg cc elementary scattering largely dominates the process up to a; ~ 0.6 (see Fig. 1), any sizeable single 
spin asymmetry measured in p^p —> DX at moderate xp 's would be the direct consequence of a non zero contribution 
of A^/g/pT • For Xp > 0.6 the competing qq — > cc term becomes approximately the same size as gg cc (see Fig.^; 
consequently the quark and gluon Sivers functions could contribute to An in approximately equal measure making 
the data analysis more involved, as we shall discuss below. 

Since we have no information about the gluon Sivers function from other experiments, we are unable to give 
predictions for the size of the An one can expect to measure at RHIC. Instead, we show what asymmetry one can 
find in two opposite extreme scenarios: the first being the case in which the gluon Sivers function is set to zero, 
A-'^ fg/pX {xa, k±a) = 0, and the quark Sivers function f^/pr {xa, k±a) is taken to be at its maximum allowed value 
at any Xa', the second given by the opposite situation, where A^/^/pt = and A^/g/p? is maximized in Xa- 

Concerning the k± dependence of the unpolarized pdf 's and the Sivers functions, we adopt, both for quarks and 
gluons, a most natural and simple factorized Gaussian parameterization 

/(x,fc±) = /(a;)-^e-'=i/<'=i>, (10) 
A-fix,k^) = A"/(x) ^e-i/C^i) cos(0.J , (11) 

where M = y/ {k'^) and (f>k^^ is the k± azimuthal angle. The extra factor 2kj_M / {k'^ + A-P) in the Sivers function is 
chosen in such a way that, while ensuring the correct small k± behaviour, it equals 1 at fcj^ = M, being always smaller 
at other values. The azimuthal cos(0/c^) dependence is the only one allowed by Lorentz invariance, via the mixed 
vector product P ■ (p x k±) where, with our frame choice, P = (0, 1, 0) is the proton polarization vector, p = (0, 0, 1) 
is the unit vector along the polarized proton motion and fcx = (cos(0fc^ ), sin(0fc^ ), 0). 
The Sivers functions (|ll|l for both quarks and gluons must respect the positivity bound 

|A^/,/pT(a;a,fej.a)| . , vy i. noA 

< 1 VXa, kj_a , (12) 



2 fa/pi^a, k±a 

which means that Eq. (|12|l can be satisfied for any Xa and k±a values by taking 

A''fa/priXa)<2fa/p{Xa). (13) 
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FIG. 2: Maximized values of |^jv| for the process p^p DX as a function of Ed and pr at fixed pseudo-rapidity (a), and as 
a function of a;^ at fixed transverse momentum (b), calculated using saturated Sivers functions, according to Eq. 1131 of the 
text. The dashed line corresponds to a maximized quark Sivers function (with the gluon Sivers function set to zero), while the 
dotted line corresponds to a maximized gluon Sivers function (with the quark Sivers function set to zero). 



For the fragmentation function Djj^q^z, k±o) we adopt a similar model, in which we assume factorization of z and 
k^D dependences 

DD/Q{z,kj_D) ^ DD/Q{z)g{ki^D), (14) 

where Du/q{z) is the usual fragmentation function available in the literature (see for instance Ref. T^) and g{k±£>) is 
a gaussian function of |fe_L_Dp analogous to that in Eq. (|10|1 . normalized so that, for a fragmenting quark of momentum 

J (fko 5{kD ■ Pc) Dd/q{z, ko) = Do/q{z) . (15) 

In Fig. ^a) we show the unpolarized cross section for the process pp DX at ^/s = 200 GeV as a function of 
both the heavy meson energy Ed and its transverse momentum pT, at fixed pseudo-rapidity 77 = 3.8 (notice that 
xp — Ed/{100 GeV)). In Fig. ^b) the same total cross section is presented as a function of xp at fixed pT = 1.5 
GeV/c. The x and Q^-dependent parton distribution functions fq/p{x, Q^) are taken from MRSTOl while the kj_ 
dependence is fixed by Eq. H1U|I with ^ {k'j_) — 0.8 GeV/c similarly, the fragmentation functions Dd/q{z, Q^) are 
from Ref. with the k± dependence fixed by ^ {k\) — 0.8 GeV/c. We have explicitely checked that our numerical 
results have very little dependence on the (fc^.) value of the fragmentation functions. Finally, we have taken as QCD 
scale = uIq. The dashed and dotted lines correspond to the qq cc and gg — > cc contributions respectively, 
whereas the solid line gives the full unpolarized cross section. These plots clearly show the striking dominance of the 
gg —^ cc channel over most of the Ejj and xp ranges covered by RHIC kinematics. 

Fig. 121 shows our estimates for the maximum value of the single spin asymmetry in p^p DX. The dashed line 
shows \An\ when the quark Sivers function is set to its maximum, i.e. f^j/pt (x) = 2fg^p{x), while setting the 
gluon Sivers function to zero. Clearly, the quark contribution to An is very small over most of the kinematic region, 
at both fixed pseudo-rapidity and varying Ejj^ Fig. 2(a), and fixed px and varying xp, Fig. 2(b). The dotted 
line corresponds to the SSA one finds in the opposite situation, when A^/^/pi (x) — 2fg/p{x) and f^^pr — 0: 
in this case the asymmetry presents a sizeable maximum in the central Ed and positive xp energy region (in our 
configuration positive xp means D mesons produced along the polarized proton direction, i.e. the positive Z-axis). 
This particular shape is given by the azimuthal dependence of the numerator of Ajv, see Eqs. ^ and 1)11(1 . When 
the energy Ed is small, pr is also very small (for instance, for Ed < 23 GeV, pt < 1 GeV/c) and the partonic 
cross sections da/dt depend only very weakly on 4>k±a- Therefore, when we integrate over (j^k^^ the partonic cross 
sections multiplied by the factor cos((/)fcj^^) from the Sivers function, we obtain negligible values. The transverse 
momentum pT of the detected D meson grows with increasing Ed and the partonic cross sections become more and 
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more sensitively dependent on (l>k±a.'- then An grows and a peak develops in correspondence of (f>k±a, — 0. Similarly, 
one can understand the behaviour of |^Ar|max in the negative and positive xp regions, Fig. 2(b). Only at very large 
Ed and xp the qq cc contribution becomes important, and a rigorous analysis in that region will only be possible 
when data from independent sources will provide enough information to be able to separate the two contributions. 

By looking at Fig.|2]it is natural to conclude that any sizeable single transverse spin asymmetry measured by STAR 
or PHENIX experiments at RHIC in the region Ed < 60 GeV or —0.2 < xp < 0.6, would give direct information on 
the size and importance of the gluon Sivers function. 



III. COMMENTS AND CONCLUSIONS 



We have shown that the observation of the transverse single spin asymmetry Aj^ for D mesons generated in p^p 
scattering offers a great chance to study the Sivers distribution functions. This channel allows a direct, uncontaminated 
access to this function since the underlying elementary processes guarantee the absence of any polarization in the final 
partonic state; consequently, contributions from Collins-like terms cannot be present to influence the measurement. 
Moreover, the large dominance of the gg — *■ cc process at low and intermediate xp offers a unique opportunity to 
measure the gluon Sivers distribution function A^/g/pt- 

Once more, intrinsic parton motions play a crucial role and have to be properly taken into account. Adopting a 
simple model to parameterize the k± dependence we have given some estimates of the unpolarized cross section for D 
meson production, together with some upper estimate of the SSA in the two opposite scenarios in which either A^/^/pT 

is maximal and A^/^/pt = or A^/^/pt = and A^/^/pt is maximal. Our results hold for D = , , D'^ , . 
Both the cross section and Ajsf could soon be measured at RHIC. 

It clearly turns out that any sizeable contributions to the p^p DX single spin asymmetry at low to intermediate 
Ed^s or xp^s would be a most direct indication of a non zero gluon Sivers function. 
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